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This  Fortran  program  generates  continuous  time  sequences  of  ground 

clutter  that  would  be  received  on  an  airborne  radar  platform  for  a  time- 

The  method  is  based  on  that  in  Reference  1. 
[2,3] 


versions 


varying  engagement  geometry. x 
This  program  replaces  e^xli^r 

Signals  are  generated  for  the  three  monopulse  channels.  It  is  assumed 
that  these  signals  will  be  radiated  into  specific  positions  on  the  receive 
antenna  pattern  so  that  each  signal  will  be  received  by  the  proper  channel 
and  rejected  by  the  others. 

This  program  was  written  for  a  single  range  gate  of  a  pulse-Doppler 
radar,  although  it  is  easily  extended  to  multiple  range  gates  (but  not 
in  real  time  with  a  single  AP120B) .  The  organization  of  this  program  is 
sketched,  iu  Figcrte.  !£x  At  the  RFSS  it  is  assumed  that  the  real-time  state 
of  the  environment  will  be  set  up  in  C0MM0N/C1/  and  periodic  calls  will 
be  made  to  CLUTTR  in  the  Datacraft/6.  The  subroutines  that  actually 
generate  the  clutter  sequences  (FILL,  CURVE,  PARAB,  CONVRT,  TIMESQ,  GAUSS4, 
FOURT)  will  be  implemented  in  the  AP120B.  Their  Fortran  equivalents  are 
included  in  this  package. 

A  test  program  (TEST)  is  also  included  to  compare  results  on  different 
computers.  The  output  of  this  program  is  a  power  spectral  estimate  of  the 
clutter  sequences  generated  under  stationary  conditions.  ^ 

The  geometry  is  shown  in  Figure  2,  with  all  angles  being  used  consistently 
throughout  the  program. 
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PROGRAM  SUMMARY 


Main  or  Driver 


TEST  (for  this  example) 


Initialization 


RCDATA 

TIMSQI 

GAUSSI  (GAUSS 3, RAND) 
Real-time  Fortran 


CLUTTR (GAIN) 
PARMS(AINT) 


Utility 


PR 

ACCUM 

POWSPT 

XMIT 

FOURT 

Real-time  AP120B 


FILL 

CURVE 

PARAB 

CONVRT 

TIMESQ 

GAUSS4 
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PROGRAM  TEST ( INPUT,  OUTPUT,  TAPE5-INPUT.  TAPE6-0UTPUT ) 

TEST  OF  SUBROUTINE  CLUTTR 

THE  TEST  CASE  CORRESPONDS  TO  ALPHA-20.  , BETA-O.  ,  GAMMA-10.  » DELTA-O.  ,  AND 
OMEGA-— 22.  5  (ALL  DEGREES). 

D I MENS ION  SR ( 2240,  3 ) ,  S I ( 2240,  3 ) ,  S ( 200,  3 ) ,  WORK ( 600 ) 

COMMON  /Cl/  WL,  H.  V,  R,  FR,  DE,  OM,  THO,  PHO 
COMMON  /C2/  PTGDSG, SIGO.  DR 
COMMON  /C3/  NFFT,  NOVLP,  NF,  A (64) 

COMMON  /C5/  PEFF, 

1  XR1  (320),  XU  (320),  XR2(320>,  XI2(320),  XR3(320),  XI3(320> 

THE  FOLLOWING  PARAMETERS  MUST  BE  SET . 

DATA  WL,  H,  R,  FR.  V/.  02.  1736.  .  10000.  ,  10000.  ,  500.  / 

DATA  DE.  OM,  THO,  PHO/O.  ,  -.  3927,  .  3244,  .1313/ 

DATA  PTGDSQ,  DR,  SIGO/1.  ,  20.  .  .  01/ 

DATA  NFFT.  NOVLP /256.  32/ 

HERE  WE  INITIALIZE 

CALL  RCDATA 
CALL  TIMSQI 
CALL  OAUSSI 

HERE  WE  GENERATE  CONTIGUOUS  TIME  SEQUENCES  IN  THE  SR  AND  SI  ARRAYS 


L-l 

DO  20  I-l. 10 
CALL  CLUTTR 

CALL  XMIT(NF,  XR1,  SR(L,  1 ) ) 

CALL  XMIT (NF, XR2.  SR (L,  2) ) 

CALL  XMIT (NF, XR3,  SR(L,  3) ) 

CALL  XMIT(NF,  XU,  SI(L,  1)  ) 

CALL  XMIT (NF,  XI2,  SI (L,  2) ) 

CALL  XMIT (NF,  XI3,  SI (L,  3) ) 

L-L+NF 
20  CONTINUE 

HERE  WE  COMPUTE  THE  (AVERAGED)  POWER  SPECTRUM 


L-l 

DO  30  1-1, 11 

CALL  PQWSPT(SR(L»  1 ).  SI  (L,  1  ) ,  200,  200,  O.  .WORK) 
CALL  POWSPT(SR(L,  2).  SI  (L,  2),  200,200,0.  .WORK) 
CALL  POWSPT(SR(L,  3), SUL,  3),  200,  200,  0.  .WORK) 
CALL  ACCUM(8R(L,  1 ) ,  S( 1,  1 ) ,  200,  I ) 

CALL  ACCUM(SR(L.  2),  S(l,  2),  200,  I) 


CALL  ACCUM<SR<L<  3) ,  S< 1 , 3) ,  200, 1 ) 
L-L+200 
30  CONTINUE 

CALL  PR  <8.200,  3) 


STOP 

END 
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SUBROUTINE  RCDATA 
C 

C  THIS  SUBROUTINE  READS  THE  CLUTTER  DATA  BASE  INTO  COMMON  /C6A 
C 

DIMENSION  IW(4),  W<22> 

COMMON  /C6/  S( 180). IS(4), A1 (216).  IA1 (36),  A2( 180),  IA2(36),  E(216), 
1  IE (36) 

DATA  IW/10HIS . »  10HIA1 . .  10HIA2 . .  10HIE . / 


DATA 

W/10HS1 . 

....  10HS2.  .  .  . 

...» 10HS3 . 

...» 10HS4 . 

1 

10HS5 . 

...»  1 OH A 1 1 .  . . 

...» 10HA12 _ 

.... 10HA13 . 

2 

10HA14.  .  .  . 

....  10HA15.  .  . 

.... 10HA1P _ 

.... 10HA22 . 

3 

10HA23.  .  .  . 

....  10HA24. . . 

....  10HA25.  .  .  . 

...» 10HA2P . 

4 

10HE1 . 

.... 10HE2 _ 

.  .  .  , 10HE3 . 

.... 10HE4 . 

5 

10HE5 . 

....  10HEP.  .  .  . 

.  .  .  / 

READ  102,  S 
READ  101.  IS 
READ  102.  A1 
READ  101,  IA1 
READ  102,  A2 
READ  101, IA2 
READ  102.  E 
READ  101, IE 

PRINT  200 

Kl=l 

K2-36 

DO  10  L»l,5 

PRINT  202,  W(L),  <S(K),  K«K1,  K2) 
Kl-Kl+36 
K2-K2+36 
10  CONTINUE 

PRINT  201. IW(1), IS 

Kl-1 

K2-36 

DO  20  L-6. 1 1 

PRINT  202, W(L).  (A1(K>.K»K1, K2) 
Kl-Kl+36 
K2-K2+36 
20  CONTINUE 

PRINT  201,  IW(2) ,  IA1 

Kl-1 

K2-36 

DO  30  L-12, 16 

PRINT  202.  W(L),  (A2(K),  K-Kl,  K2) 


Kl-Kl+36 
K2-K2+36 
30  CONTINUE 

PRINT  201,  IW<3>,  IA2 

Kl=»l 

DO  40  L*17,22 

PRINT  202,  W<L),  <E<K).  K=Kl,  K2) 

Kl»Kl+36 
K2-K2+36 
40  CONTINUE 

PRINT  201,  IW(4) » IE 

RETURN 

100  FORMAT < 20F4.  0) 

101  FORMAT <3612) 

102  FORMAT < 2< 9F4.  0,  4X) ) 

200  FORMAT < 1H1. 50HGR0UND  CLUTTER  DATA  BASE  READ  BY  SUBROUTINE  RCDATA/ ) 

201  F0RMAT</1XA1 0,3X3613) 

202  FORMAT < /I XA10.  3X9F7.  3/U4X9F7.  3)) 

END 
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subroutine  timsgi 

THIS  SUBROUTINE  PERFORMS  INITIALIZATION  FOR  SUBROUTINE  TIMESG. 

THE  INITIALIZATION  PARAMETERS  ARE . 

NFFT  -  SIZE  OF  FFT  < «2*# I NTEGER > 

NOVLP  *  NUMBER  OF  SAMPLES  THAT  THE  FFTS  OVERLAP  < NFFT /B  IS  OK> 

THE  ARRAYS  MUST  BE  DIMENSIONED  AS  LARGE  AS . 

XR 1 . X 1 1 .  ETC  NFFT +NOVLP-1 

A  NOVLP -1 

COMMON  /C3/  NFFT,  NOVLP,  NF,  A (64 > 

COMMON  /C5/  PEFF, 

1  XR1<320),  XII  <320),  XR2<320),  XI2(320>,  XR3<320),  XI3O20) 

NF-NFFT-NOVLP 
Nl-NOVLP-1 
N2-NFFT+N1 
DO  lO  K-1.N1 
A  <  K  > -K /FLOAT  <  NOVLP ) 

10  CONTINUE 

CALL  XMIT (— N2,  0.  ,  XR1 ) 

CALL  XMIT <  — N2.  O.  .  XR2) 

CALL  XMIT < -N2. O.  .  XR3) 

CALL  XMIT<-N3, O.  ,  XII) 

CALL  XMIT (— N2,  0.  ,  XI2) 

CALL  XMIT < -N2,  0.  ,  XI3) 

RETURN 

END 
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SUBROUTINE  GAUSS I 


THIS  SUBROUTINE  FILLS  ARRAYS  TR  AND  TI  WITH  INDEPENDENT  GAUSSIAN 
RANDOM  NUMBERS  FOR  USE  WITH  SUBROUTINE  GAUSS4. 

ARRAYS  TR  AND  TI  MUST  BE  DIMENSIONED  AS  LARGE  AS  NTR+NMAX  AND  NTI+NMAX 
RESPECTIVELY. 

IF  N  BELOW  IS  NOT  EVEN,  ADD  1  TO  THE  DIMENSION  OF  TI. 

COMMON  /CS/  IRSET,  NTR,  NTI, NMAX,  KR,  KI,  TR C 1253),  TI< 1353) 

DATA  NTR,  NTI/997,  1097/,  NMAX/256/,  KR.KI/1,1/ 

N*NTR+NTI+2*NMAX 
DO  20  K=1,N,  2 

CALL  GAUSS3 i TR ( K ) ,  TR ( K+ 1 ) ,  1 .  ) 

20  CONTINUE 
RETURN 
END 
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SUBROUTINE  0AUSS3( X, Yi  P ) 

THIS  SUBROUTINE  GENERATES  RANDOM  PAIRS  FOR  THESE  DISTRIBUTIONS . 

CALL  0AUSS3<X,  Y,  P) - GAUSSIAN  COMPONENTS  OF  AVERAGE  POWER  P 

CALL  RANG3(  X,  Y»  P ) . RANDOM  PHASOR  COMPONENTS  OF  POWER  P 

DATA  I/O.  / 

C  ) X ( FNAR® ) X  <  DNAR 

IF(P )  30,20,10 
10  E*SGRT  < -P*ALOG  <  RAND ( Z ) ) ) 

00  TO  15 
ENTRY  RAN03 
IF(P )  30,20,12 
12  E-SQRT<P> 

19  A-RAND(Z) 

A»A+A-1. 

B-RAND(Z) 

B-B+B-l. 

A2=A*A 

B2*B*B 

C-A2+B2 

IF (C.  GT.  1.  )  GO  TO  15 
X=E*<A2-B2)/C 
Y-E*2.  *A*B/C 
RETURN 

20  X-0. 

Y*0. 

RETURN 

30  STOP  44 
END 
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FUNCTION  RAND ( Z ) 

C  FOR  TEST  PURPOSES  ONLY. 

DIMENSION  A( 8) 

DATA  A/.  1»  .  2»  .  3»  .  4.  .  5»  .  .  7.  .  8/ 

DATA  K/l/ 

R AND= A ( K ) 

K-K+l 

IF(K.  CT.  8)  K=1 

RETURN 

END 
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SUBROUTINE  CLUTTR 

IN  THIS  SUBROUTINE  WE  GENERATE  THE  CLUTTER  MODULATION  SIGNALS  FOR  ONE 
RANGE  GATE  ON  EACH  OF  THREE  MONOPULSE  CHANNELS. 

THE  BACKSCATTER  COEFFICIENT  OF  THE  GROUND  IS  ASSUMED  TO  VARY  AS  THE 
SINE  OF  THE  GRAZING  ANGLE. 

THE  REAL-TIME  PARAMETERS  ON  INPUT  ARE . 

WL  -  WAVELENGTH 
H  -  PLATFORM  ALTITUDE 
V  -  PLATFORM  VELOCITY 

R  -  RANGE  TO  RANGE  GATE  OF  INTEREST  (UNAMBIGUOUS) 

FR  =  EQUIVALENT  SAMPLE  RATE  IN  RECEIVER  (PROCESSING  BAND) 

DE  -  DELTA,  THE  MISSILE  DIVE  ANGLE 

OM  -  OMEGA,  THE  MISSILE  ROLL  ANGLE  (CCW  IS  POSITIVE) 

THO  -  THETA-ZERO,  THE  AZIMUTH  ANGLE  OF  THE  ANTENNA  BORESIGHT 
PHO  -  PHI-ZERO,  THE  ELEVATION  ANGLE  OF  THE  ANTENNA  BORESIGHT 

THE  NONREAL-TIME  PARAMETERS  ON  INPUT  ARE . 

PTGDSQ  =  PRODUCT  OF  TRANSMIT  POWER,  GAIN,  SQUARE  OF  CHAMBER  LENGTH 
DR  -  RANGE  RESOLUTION 

SIGO  -  GROUND  BACKSCATTER  COEFFICIENT  AT  30-DEG  GRAZING  ANGLE 

THE  REAL-TIME  OUTPUT  FOR  THAT  PORTION  COMPUTED  IN  THE  DATACRAFT  IS 
THROUGH  /C4/»  AND  IT  IS  COMPUTED  IN  SUBROUTINE  PARMS. 

SOME  OF  THE  OTHER  QUANTITIES  USED  IN  THIS  SUBROUTINE  ARE . 

AL  =  ALPHA,  THE  AZIMUTH  ANGLE  OF  THE  BEAM  AXIS  IN  GROUND  COORD. 
BE  -  BETA,  THE  ELEVATION  ANGLE  OF  THE  BEAM  AXIS  IN  GROUND  COORD. 
GA  -  GAMMA,  THE  GRAZING  ANGLE  AT  THE  RANGE  RING  ON  THE  GROUND 

ALL  ANGLES  ARE  IN  RADIANS  AND  ALL  DISTANCES  IN  METERS. 

COMMON  /Cl/  WL,  H,  V,  R, FR, DE, OM, THO,  PHO 
COMMON  /C2 /  PTGDSQ,  SIGO, DR 
COMMON  /C3/  NFFT 
DATA  FCTR/1.  / 

COMPUTE  GEOMETRY  AND  DOPPLER  PARAMETERS 
SGA-H/R 

CGA-SQRTd.  -SQA**2> 

C A- AT AN < SC A /CO A) 

CPHO-COS(PHO) 

SPHO»SIN(PHO> 

CTHO=C  OS ( THO ) 
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STHO—S I N  <  THO ) 

COM-COS(OM) 

SOM-SIN(OM) 

CDE-COS(DE) 

SDE—SIN1DE) 

AX-CPHO*STHO 

AY«CPHO*CTHO 

AZ-SPHO 

BX«AX*COM-AZ*SOM 

BY-AY 

BZ*AX*SOM+AZ*COM 
CX»  BX 

CY»  BY*CDE+BZ*SDE 
CZ«-BY*SDE+BZ*CDE 
AL-ATAN<CX/CY> 

BE-ATAN<-CZ/SQRT<1.  -CZ**2>> 

PEFF-.  0 1 3*PT GDSQ* < 2.  *SICO*SGA ) *DR*GAIN< GA-BE) /R**3 
FMAX-2.  *V/WL 

FMC  -FMAX*  <  CCA*CDE+SGA*SDE ) 

FSCL-.  9848*FMAX/FMC 
IAMB—FMC/FR 

FNP1— ( IAMB-1 )*FR*FSCL/FMAX 
DFNP-FNP1 /FLOAT <  < IAMB-1 )*NFFT> 

FNP1=.  9848-FCTR*(.  9848-FNP1) 

DFNP-FCTR*DFNP 
ALD=57.  3*AL 
OMD-57.  3*0M 
IF<  ALD.  OT.  0.  )  GO  TO  10 
0MD-360.  -OMD 
ALD—ALD 
lO  OMDP-OMD 

IFIOMD.  GT.  315.  )  DMDP-OMDP-360. 

ICASE—  ( OMDP+45.  )/90. 

OMDP-OMDP - 1 C ASE *90 . 

I ROT-O 

IF(  I  CASE.  EG.  1.  OR.  ICASE.  EG.  3)  IROT-1 

)/4.  21E/  4.  9F3,  1.  9F2,  1.  BF3,  4.  8F3,  413/ . RTTULC  H21/ CTAMROF  001 

FFEP,  PNFD,  1PNF  1 

. LCSF. CMF,  XAMF.  PDMO, DMO, DLA, EB. LA, AG, TORI, ESACI,  BMAI,  001  TNIRP 
COMPUTE  REAL-TIME  PARAMETERS 

CALL  P ARMS  <  ALD,  OMDP ,  FNP 1 ,  DFNP ) 

GENERATE  REAL-TIME  SPECTRA 
IROT-O 

IF<  ICASE.  EG.  1.  OR.  ICASE.  EG.  3)  IROT-1 
CALL  FILLUROT) 

RETURN 

END 


m 
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FUNCTION  GAIN (PHI) 

HERE  WE  COMPUTE  THE  2-WAY  POWER  GAIN  BY  TABLE  LOOKUP.  PHI  IS  THE 
ANOLE  IN  RADIANS. 

NB  =  NUMBER  OF  SAMPLES  IN  TABLE 
DB  «  SAMPLE  SPACING  (DEG) 

THE  FIRST  SAMPLE  IS  AT  PH 1=0. 

DIMENSION  BEAM(lOl) 

DATA  NB/101/j  DB/.  5/ 

DATA  A/.  00389/.  B/2.  30/ 

DATA  NO/O/ 

IF(NO.  GT.  O)  GO  TO  30 
N0=1 

DO  20  K=1,NB 

BEAM ( K ) =EXP ( -A.  *A* ( ( K- 1 ) *DB >  **B ) 

20  CONTINUE 
30  E=S7.  3*ABS(PHI) 

IE-E/DB+1.  5 
GAIN-0. 

IF(IE.  LE.  NB)  GAIN=BEAM( IE) 

RETURN 

END 
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SUBROUTINE  PARMS< ALPHA,  OMEGA,  FI. DF) 

IN  THIS  SUBROUTINE  WE  COMPUTE  THE  REAL-TIME  PARAMETERS  THAT  WILL  BE 
SENT  TO  THE  AP120B.  THE  PARAMETERS  ARE  DETERMINED  BY  INTERPOLATION 
WITHIN  THE  DATA  BASE  STORED  IN  COMMON  /C6/. 

THE  INPUT  IS . 

ALPHA  -  ABSOLUTE  VALUE  OF  AZIMUTH  ANGLE  (DEG) 

OMEGA  -  MISSILE  ROLL  ANGLE  (DEG,  -45.  .  LE.  OMEGA.  LE.  45.  ) 

FI  -  NORMALIZED  DOPPLER  OF  FIRST  SAMPLE 
DF  -  SAMPLE  SPACING  IN  DOPPLER 

THE  OUTPUT  IS  THROUGH  COMMON  /C4/ . 


IS,  SPK . SUM  SPECTRUM  PARAMETERS 

IA1.A1PK _ DEL-AZ  SPECTRUM  PARAMETERS,  LEFT  HALF 

IA2, A2PK _ DEL-AZ  SPECTRUM  PARAMETERS,  RIGHT  HALF 

IE. EPK .  DEL-EL  SPECTRUM  PARAMETERS 

PSL . SIDELOBE  LEVEL  (SEE  SUBROUTINE  FILL) 


SKP, A1PK, A2PK,  EPK,  AND  PSL  MUST  BE  IN  THE  SAME  UNITS  (NEPERS). 

TNI A  LANRETXE 

COMMON  /C4/  IS(6),  IA1 (6),  IA2(6>,  IE(6>, SPK,  A1PK,  A2PK,  EPK,  PSL 
COMMON  /C6/  SI (36) ,  S2(36) ,  S3(36) , S4(36) ,  55(36) ,  IST(4) , 

1  All (36).  A12 ( 36 ) ,  A13( 36) ,  A14(36) , A15(36 ) ,  A1P ( 36) ,  I AIT (36), 

2  A22 ( 36 ) ,  A23 ( 36 ) ,  A24 ( 36 ) ,  A25 ( 36 ) ,  A2P ( 36 ) ,  IA2T(36), 

3  E 1 ( 36 ) ,  E2 ( 36 ) ,  E3 ( 36 ) , E4 ( 36 ) , E5 ( 36 ) ,  EP ( 36 ) ,  IET(36> 

COMMON  /CINT/  Kl,  K2,  K3,  K4,  HOM,  HAL 

DATA  NAL.  NOM/4,  9/ 

DATA  DAL,  DOM/ 10.  ,11.  25/ 

ISAMP(A)«(AINT(A)-F1)/DF+1.  5 

HERE  WE  FIND  THE  PARAMETERS  FOR  INTERPOLATION 

AL-ALPHA/DAL 

IAL-AL 

HAL-AL-IAL 

IAL-IAL+1 

IF(IAL.  OE.  NAL)  HAL-1. 

IFdAL.  OE.  NAL)  IAL-NAL-1 

OM- ( OMECA+45.  )/DOM 

IOM-OM 

HOM-OM-IOM 

IOM-IOM+1 

IF(IOM.  GE.  NOM)  HOM-1. 

IF ( IOM.  GE.  NOM)  IOM-NOM-1 
Kl-( IAL-1 )*NOM+IOM 
K2-K1+1 
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K3-K1+N0M 
K4-K2+N0M 
KAL-ALPHA/DAL+1.  5 
KAL«MI NO  (  KAL# NAL ) 

KOM=  (  OMEGA+4  5.  )/DOM+l.  5 
K»(KAL-1)*N0M+K0M 

HERE  WE  COMPUTE  THE  REAL-TIME  PARAMETERS 

IS(1)  =*ISAMP  (SI ) 

IS(2)  «ISAMP(S2> 

ISO)  *ISAMP(S3) 

IS(4)  -ISAMP(S4) 

18(5)  -ISAMP(S5) 

IA1 < 1 )*ISAMP ( A1 1 ) 

IA1 (2)*ISAMP(A12> 

IA1(3)*ISAMP(A13) 

IA1(4)«ISAMP(A14) 

IA1 <5>»ISAMP ( A15) 

IA2(1)-IA1(5) 

IA2(2)*ISAMP ( A22) 

IA2(3)«ISAMP(A23> 

I A2 ( 4 ) * I SAMP  <  A24 ) 

IA2<  5)“ISAMP  < A25) 

IE( 1 )  -ISAMP(El) 

IE(2)  ®ISAMP(E2) 

IE(3)  -ISAMP(E3) 

IE(4)  «ISAMP(E4) 

IE( 5)  *ISAMP(E5) 

SPK  -40. 

PSL  -0. 

A1PK*AINT(A1P) 

A2PK*A I NT ( A2P  > 

EPK  «AINT(EP) 

18(6)  -IST(KAL) 

IA1(6)«IA1T(K> 

IA2(6)-IA2T(K> 

IE(6)  «IET(K) 

CONVERT  DB  TO  NEPERS  (AMPLITUDE) 

SPK  11513*SPK 
A1PK-.  1 1513*A1PK 
A2PK*.  11513*A2PK 
EPK  «.  11513*EPK 
PSL  «.  11513*PSL 


RETURN 

END 
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FUNCTION  AINT< A) 

INTERPOLATION  FOR  SUBROUTINE  PARMS 
DIMENSION  A( 1 > 

COMMON  /CINT/  Kl»  K2»  K3»  K4*  HOMi  HAL 
Al*(l.  -H0M)*A(K1)+H0M*A(K2) 

A2*  < 1 .  -HOM  >  *A  <  K3 ) +HOM*A ( K4  > 

A I NT*  < 1 .  -HAL ) *A 1 +HAL*A2 

RETURN 

END 


J 
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SUBROUTINE  PR  (A,  NX.  NY) 

IN  THIS  SUBROUTINE  WE  PRINT 
EQUIVALENT  TO  THE  DIMENSION 

THE  PRINTING  IS  ON  FILE  LU, 


THE  2-DIMENSIONAL  ARRAY  A  WITH  THE  FORMAT 

A (NX, NY> . 

WHICH  MUST  BE  EQUIPPED. 


DIMENSION  A( 1 > 

COMMON  /MESGE/  MFLAG 
DATA  LU/6/ 

IF ( MFLAG.  NE.  1)  WRITE  (LU.  100) 

MFLAG=0 

DO  20  J-l.NY 

L1*(J-1)*NX+1 

L2-L1+9 

L2MAX-J*NX 

L2»MIN0(L2.  L2MAX) 

DO  20  1*1,  NX, 10 

IFd.LE.  1)  WRITE  (LU.  101)  J,  (A(K> , K*L1,  L2> 
IF(I.CT.  1)  WRITE  (LU.  102)  ( A(K> .  K=L1,  L2> 
Ll-Ll+10 

L2*MIN0(L2+10,  L2MAX) 

20  CONTINUE 
RETURN 

100  FORMAT ( 1H1 ) 

101  F0RMAT(/I4,  2X10E12.  4) 

102  FORMAT (6X10E 12.  4) 

END 


i 
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SUBROUTINE  ACCUM< X,  Y,  N,  1REP ) 


IN  THIS  SUBROUTINE  WE  ACCUMULATE  AN  ARRAY. 
REPRESENTS  A  RUNNING  AVERAGE. 

DIMENSION  X(1)»Y(1) 

IF <  IREP.  EQ.  1)  GO  TO  30 
REP- I REP 
DO  20  K*l*  N 

Y(K>  — <  (REP— 1.  >*Y(K>+X(K> )/REP 
20  CONTINUE 
RETURN 

30  CALL  XMIT <N#  X»  Y> 

RETURN 

END 


AT  ALL  TIMES  ARRAY  Y 
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SUBROUTINE  POWSPT<XR, XI.  NIN,  NOUT,  ALPHA.  W) 

IN  THIS  SUBROUTINE  WE  COMPUTE  THE  POWER  SPECTRUM  OF  THE  COMPLEX  TIME 
SEQUENCE  IN  THE  ARRAY-PAIR  <XR,XI>  OF  LENGTH  NIN.  THE  POWER  SPECTRUM 
IS  RETURNED  IN  ARRAY  XR.  AND  IT  IS  NOW  OF  LENGTH  NOUT. 

THE  SAMPLE  SPACING  OF  THE  POWER  SPECTRUM  IS  1/NOUT  OF  THE  REPETITION 
FREQUENCY. 

A  COSINE-ON-A-PEDESTAL  WEIGHTING  IS  APPLIED  TO  THE  INPUT  SAMPLES. 
ALPHA  IS  THE  RATIO  OF  THE  WEIGHTING  FUNCTION  AT  THE  EDGE  TO  THE 
CENTER.  ALPHA-.  08  FOR  HAMMING  AND  ALPHA- 1.  0  FOR  UNIFORM  WEIGHTING. 

ARRAY  W  IS  A  WORKING  ARRAY  AND  IT  MUST  BE  DIMENSIONED  AS  LARGE  AS . 

NIN  IF  NOUT.  EQ.  2** INTEGER 

NIN+2*N0UT  IF  NOUT.  NE.  2** INTEGER 

THE  WEIGHTS  ARE  NORMALIZED  SO  THAT  THE  SUM  IS  UNITY. 

DIMENSION  XR(1).XI(1).W(1) 

DATA  TWOPI/6.  2831853/ 

DATA  NO/O/ 

IF  (NO.  GT.  0)  GO  TO  25 
NO-1 

N-NOUT-NIN 
A— ( 1.  +ALPHA) /2. 

B-<1.  -ALPHA )/2. 

CN— (NIN+l )/2. 

XN— NIN 

DO  10  K-l.NIN 

W  <  K ) =A+B*COS  <  TWOP I *  <  K-CN ) / XN  > 

10  CONTINUE 
SUM-0. 

DO  15  K-l.NIN 
SUM-SUM+W  <  K ) 

15  CONTINUE 
WNORM— SUM 
DO  20  K-l.NIN 
W  <  K ) -W  <  K  > /WNORM 
20  CONTINUE 
25  IF(N)  60.35.30 
30  CALL  XMIT (-N.  0.  . XR< NIN+l ) > 

CALL  XMIT(-N.  0.  .  XI (NIN+l ) ) 

35  DO  40  K-l.NIN 
XR(K)-XR(K)*W(K) 

XI<K)-XI(K)*W<K) 

40  CONTINUE 

CALL  FOURT ( XR.  XI,  NOUT,  1,  1, 1. W(NIN+1 >» W(NIN+N0UT+1 > > 

DO  50  K-l.NOUT 
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XR<K)«XR(K)**2+XI<K>**2 
50  CONTINUE 
RETURN 

60  PRINT  100, NIN.NOUT 
STOP 

100  FORMAT  < / 1 X49H***ERRGR  IN  PGWSPT  .  NOUT  IS  SMALLER  THAN  NIN  *##/ 
t  /10X4HNIN*I6,  10X5HN0UT*I6> 

END 


<J<J<JU<J<JU 


SUBROUTINE  XKIT (N»  A. B ) 

IN  THIS  SUBROUTINE  WE  EITHER  TRANSMIT  ARRAY  A  TO  ARRAY  B  (IF  N. GT.  0> 
OR  WE  TRANSMIT  THE  CONSTANT  A  TO  ARRAY  B  (IF  N.  LT.  0>.  IN  EITHER  CASE 
TIJF  ARRAY  LENGTH  IS  IABS(N). 

THIS  SUBROUTINE  SHOULD  BE  WRITTEN  IN  ASSEMBLY  LANGUAGE 

DIMENSION  A(1>.B(1) 

IF(N>  10/20,25 
10  NN=-N 
AA— A ( 1 ) 

DO  15  K-l.NN 
B(K>-AA 
15  CONTINUE 
20  RETURN 
25  DO  30  K-l.N 
B(K)*A(K) 

30  CONTINUE 
RETURN 
END 


OflO 
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SUBROUTINE  FILL< IROT) 

C 

C  IN  THIS  SUBROUTINE  WE  OENERATE  THE  REAL-TIME  CLUTTER  MODULATION 
C  SIGNALS  ON  THE  THREE  MONOPULSE  CHANNELS  (FOR  ONE  RANGE  GATE  ONLY). 

C  THIS  SUBROUTINE.  AS  WELL  AS  SUBROUTINES  CURVE,  CONVRT,  PARAB,  TIMESG, 

C  GAUSS4,  AND  FFT2,  WILL  BE  IMPLEMENTED  ON  THE  AP120B. 

C 

C  THE  REAL-TIME  INPUT  IS  THROUGH  /C4/ . 

C 

C  IS,  SPK . SUM  SPECTRUM  PARAMETERS 

C  IA1.A1PK _ DEL-AZ  SPECTRUM  PARAMETERS,  LEFT  HALF 

C  IA2, A2PK _ DEL-AZ  SPECTRUM  PARAMETERS,  RIGHT  HALF 

C  IE, EPK . DEL-EL  SPECTRUM  PARAMETERS 

C  PSL . SIDELOBE  LEVEL  (SEE  BELOW) 

C 

C  THE  REAL-TIME  OUTPUT  IS  THROUGH  /C5/ . 

C 

C  PEFF  -  EFFECTIVE  POWER  RADIATED  FROM  RFSS  ARRAY  (FROM  SUB.  CLUTTR ) 
C  XR 1  *  REAL  ARRAY  CONTAINING  SUM  CHANNEL  MODULATION  SIGNAL 

C  XII  *  IMAG  ARRAY  CONTAINING  SUM  CHANNEL  MODULATION  SIGNAL 

C  XR2  =  REAL  ARRAY  CONTAINING  DAZ  CHANNEL  MODULATION  SIGNAL 

C  XI2  =  IMAG  ARRAY  CONTAINING  DAZ  CHANNEL  MODULATION  SIGNAL 

C  XR3  *  REAL  ARRAY  CONTAINING  DEL  CHANNEL  MODULATION  SIGNAL 

C  XI3  -  IMAG  ARRAY  CONTAINING  DEL  CHANNEL  MODULATION  SIGNAL 

C 

C  IROT  IS  A  ROTATION  FLAG.  IF  IROT.  GT.  0  THE  AZ  AND  EL  ARRAYS  ARE 
C  INTERCHANGED. 

C 

C  THIS  PROGRAM  GENERATES  A  CONSTANT-LEVEL  SIDELOBE  SHOULDER  OF  (SPK-PSL) 
C  DB  BELOW  THE  SUM  SPECTRUM  PEAK. 

C 

C  ARRAYS  S,  A,  AND  E  ARE  FOR  TEMPORARY  STORAGE  OF  THE  SUM,  AZ-  AND  EL- 
C  DIFFERENCE  SPECTRA.  THEY  MUST  BE  DIMENSIONED  AS  LARGE  AS  2*NFFT. 

C 

C  THE  DIMENSIONED  ARRAYS  IN  COMMON  /C5/  MUST  ACCOMMODATE  NFFT+NOVLP 
C  WORDS  (SEE  TIMESQ). 

C 

COMMON  /C3/  NFFT 

COMMON  /C4/  IS(6).  IA1 (6), IA2<6>,  IE<6>, SPK, A1PK, A2PK,  EPK,  PSL 
COMMON  /C5/  PEFF, 

1  XRK320),  XII  (320),  XR2(320>,  XI2<320>,  XR3(32 0),  XI3(320> 

COMMON  /C7/  S( 512) ,  A ( 512 ) ,  E ( 512) 

COMMON  /LIM/  NFFT2 
NFFT2=NFFT  *2 

INITIALIZE  ARRAYS  TO  SIDELOBE  LEVEL 

CALL  XMIT(-NFFT2,  PSL,  S) 

CALL  XMIT(-NFFT2,  PSL,  A) 

CALL  XMIT(-NFFT2,  PSL,  E) 


non  o  n  o  o  o  o 
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GENERATE  LOO-AMPLITUDE  SPECTRA 

CALL  CURVE (St IS, SPK) 

CALL  CURVE (At IAltAlPK) 

CALL  CURVE  (At  IA2,  A2PK) 

CALL  CURVE (Et  IEt  EPK) 

CONVERT  TO  LINEAR  AMPLITUDE  AND  FOLDOVER 
SN0RM=0. 

CALL  CONVRT ( S,  SNORM ) 

CALL  CONVRT (At SNORM) 

CALL  CONVRT (E, SNORM) 

GENERATE  TIME  SEQUENCES 

CALL  TIMESQ(S. XR1,  XII ) 

IF ( IROT.  GT.  0)  GO  TO  35 
CALL  TIMESQ( At  XR2.  XI2) 

CALL  TIMESQ(E. XR3t XI3) 

GO  TO  40 

35  CALL  TIMESQ(A.  XR3»  XI3) 

CALL  TIMESQ(Ei  XR2,  XI2) 

40  RETURN 
END 


o  o  o  non  nnnnnnononnnnrjoooonnonor} 
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SUBROUTINE  CURVE(S, II, PK) 

THIS  SUBROUTINE  GENERATES  A  CURVE  IN  ARRAY  S  BETWEEN  SAMPLES  II (1) 
11(5).  THE  CURVE  IS  DEFINED  IN  FOUR  SEGMENTS  BETWEEN  FIVE  POINTS 
WHERE  THE  POINTS  ARE  GIVEN  BY . 


II  (1) 

11(2) 

11(3) 

11(4) 

11(5) 

II(6)-ITYP  IS  A  PARAMETER  THAT 

ITYP 

0 

1 

2 


VALUE 

0. 

PK/2. 

PK 

PK/2. 

0. 

DEFINES  THE  TYPE  OF  SEGMENTS . 

SEGMENTS 

LINE, LINE,  LINE,  LINE 
LINE, PARAB,  LINE,  LINE 
LINE, PARAB,  PARAB,  LINE 


THE  PARABOLA  IS  ALWAYS  CONSTRAINED  TO  HAVE  ZERO  SLOPE  AT  THE  CENTER 
POINT. 


DIMENSION  S(l). 11(6) 
COMMON  /ABC/  A,  B.  C 


IF(PK.  LE.  0.  )  RETURN 
P2-.  5*PK 
ITYP=I 1(6) 


COMPUTE  FIRST  SEGMENT 


C*0. 

B*P2/ AMAXO (11(2)— I 1(1), I) 

A— B*II(1) 

CALL  PARAB(S, 11(1),  II (2) ) 

COMPUTE  SECOND  SEGMENT 

IF( ITYP.  OT.  0)  GO  TO  20 
C-0. 

B='P2/ AMAXO  ( II(3>  — 11(2),  1) 
A»P2-B*II(2> 

GO  TO  25 

20  C— P2/AMAX0(II(3)-II(2>,  1)**2 
B«-2.  *C*I  1(3) 
A-P2-B*II(2)-C*II(2)##2 
25  CALL  PARAB (S,  11(2)+!,  11(3)) 


AND 
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COMPUTE  THIRD  SEGMENT 

IF <  ITYP.  QT.  1)  GO  TO  30 
C*0. 

B— P2/AMAX0(II<4)-II<3),  1) 
A«PK-B*II<3> 

GO  TO  35 

30  C— P2/AMAX0<II<4)-II<3),  1)**2 
B— 2.  *C*II<3) 
A»PK-B*II(3)-C*II<3>**2 
35  CALL  PARAB(S/  1 1 (3)+l»  II(4>) 

COMPUTE  FOURTH  SEGMENT 


C-0. 

B— P2/AMAX0(II(5)-II<4).  1) 
A“P2— B*I 1(4) 

CALL  PARABtS.  II<4)-H.  11(5)) 


RETURN 

END 
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SUBROUTINE  PARAB(S/  III* 112) 

THIS  SUBROUTINE  GENERATES  A  PARABOLA  OF  THE  FORM . 

S<I>  -  A  +  B*I  +  C*I#*2 

DIMENSION  S(l) 

COMMON  /ABC/  A,  B,C 
COMMON  /LIM/  LIM 
II— MAXO(IIl« 1) 

I2=MIN0(II2.  LIM) 

IFU1.0T.  12)  RETURN 
DO  20  1*11.  12 
S(I)=(C*I+B)*I+A 
20  CONTINUE 
RETURN 
END 


SUBROUTINE  CONVRT<S,  SNORM) 

C 

C  IN  THIS  SUBROUTINE  WE  CONVERT  THE  SPECTRUM  IN  ARRAY  S  OF  LOO-AMPLITUDE 
C  UNITS  TO  LI NEAR -AMPLITUDE  UNITS.  IN  ADDITION/  TWO  PRF  INTERVALS  ARE 
C  FOLDED  (SUMMED).  THE  INPUT  S-ARRAY  IS  OF  LENGTH  2*NFFT  AND  THE  OUTPUT 
C  S-ARRAY  IS  OF  LENOTH  NFFT.  AN  APPROXIMATION  IS  USED  IN  THE  DO- 20  LOOP 
C  FOR  A  FASTER  COMPUTATION. 

C 

C  FOR  THE  SUM  SPECTRUM  SET  SNORM=ZERO  ON  INPUT  (DO  NOT  USE  THE  CONSTANT 
C  O.  »  HOWZ-VER).  THE  NORMALIZATION  FACTOR  SNORM  WILL  BE  RETURNED.  FOR 
C  THE  DIFFERENCE  SPECTRA  USE  THIS  SAME  VALUE  ON  INPUT. 

C 

DIMENSION  S(l) 

COMMON  /C3/  NFFT 
DO  20  K-l.NFFT 
S(K)«EXP(S(K)+S(K+NFFT) ) 

20  CONTINUE 

IF  (SNORM.  OT.  0.  )  GO  TO  35 
SUM-0. 

DO  30  K-l.NFFT 
SUM— SUM+S ( K  >  **2 
30  CONTINUE 

SNORM- 1.  /SORT (SUM) 

35  DO  40  K-l.NFFT 
S ( K ) — S ( K ) *SNORM 
40  CONTINUE 
RETURN 
END 


ooonoftnftnftoonnonnoortoonftonoortoftftnoo 
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SUBROUTINE  TIMESQ(SA,  XR, XI  > 

IN  THIS  SUBROUTINE  WE  GENERATE  A  CORRELATED  TIME  SEQUENCE  THAT  HAS  AN 
AMPLITUDE  SPECTRAL  FUNCTION  DESCRIBED  BY  ARRAY  SA  OF  LENGTH  NFFT.  THE 
METHOD  USED  IS  THE  ON-LINE  FFT  APPROACH  DESCRIBED  IN  SECTION  6  OF  MR I 
REPORT  132-44.  WHERE  SUCCESSIVE  CALLS  TO  TIMESQ  WILL  CREATE  CONTIGUOUS 
TIME  SEQUENCES  BASED  ON  OVERLAPPING  FFTS  IN  THE  FREQUENCY  DOMAIN. 

ON  INPUT . 

NFFT  =  SIZE  OF  FFT  < -2** INTEGER ) 

NOVLP  -  NUMBER  OF  SAMPLES  THAT  THE  FFTS  OVERLAP  (NFFT/B  IS  OK) 

SA(K)  *  SPECTRAL  AMPLITUDE  OF  KTH  SAMPLE.  K-1,...,NFFT 

ON  OUTPUT . 

XR(K)  *  KTH  SAMPLE  OF  IN-PHASE  SIGNAL,  K=1 . NF 

XI (K)  »  KTH  SAMPLE  OF  QUADRATURE  SIGNAL,  K-l, . . . , NF 

WHERE  NF-NFFT-NOVLP. 

THE  ARRAYS  MUST  BE  DIMENSIONED  AS  LARGE  AS . 


SA 

NFFT 

XR 

NFFT+NOVLP-1 

XI 

NFFT +N0VLP-1 

A 

NOVLP -1 

NOTE  THAT  THE  INFORMATION  IS  STORED  IN  THE  XR  AND  XI  ARRAYS  FROM 
NFFT+1  TO  NFFT+NOVLP+1  THAT  WILL  BE  USED  ON  THE  NEXT  CALL.  PRIOR  TO 
THE  FIRST  CALL  TO  TIMESQ  THE  INITIALIZATION  SUBROUTINE  (TIMSQI)  MUST 
BE  CALLED. 

THE  WORK  ARRAYS  WRKR  AND  WRKI  NEED  BE  DIMENSIONED  ONLY  IF  NFFT.  NE. 
2<<INTEGER<  IN  WHICH  CASE  THE  DIMENSION  SIZE  MUST  BE  NFFT. 

DIMENSION  WRKR< 160),  WRKI <160) 

DIMENSION  5A(  1),XR(1),XI(1) 

COMMON  /C3/  NFFT,  NOVLP,  NF, Atl) 

CALL  CAUSS4  <  SA,  XR,  X I ,  NFFT ) 

CALL  FOURT ( XR,  X I ,  NFFT ,  1,0,  1.  WRKR,  WRKI) 

Nl-NOVLP-1 
DO  20  K-l.Nl 

XR ( K ) -A <  K  >  *  <  XR  <  K ) -XR  <  K+NFFT ) ) +XR  <  K+NFFT ) 

XI <K)-A(K)*<  XI (K)-XI (K+NFFT) ) +XI (K+NFFT) 

XR ( K+NFFT ) -XR  <  K+NFFT-NOVLP ) 

X I <  K+NFFT ) -X I <  K+NFFT-NOVLP ) 

20  CONTINUE 
RETURN 
END 
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SUBROUTINE  GAUSS4(A.  XR»  XI#  N> 

C 

C  GENERATES  RANDOM  PHASOR  COMPONENTS  < XR <K> .  XI <K> )  OF  ZERO  MEAN  WITH  THE 
C  AVERAGE  POWER  OF  XR<K)**2+XI <K)**2  OIVEN  BY  A(K)**2  FOR  K»1,...,N. 

C 

C  THE  RANDOM  NUMBERS  ARE  ACCESSED  FROM  RECIRCULATING  TABLES.  NTR  AND 
C  NTI  ARE  THE  CIRCULATION  PARAMETERS  AND  NMAX-MAX(N).  WE  SHOULD  CHOOSE 
C  NTR  AND  NTI  RELATIVELY  PRIME. 

C 

C  ARRAYS  TR  AND  TI  MUST  BE  DIMENSIONED  AS  LARGE  AS  NTR+NMAX  AND  NTI+NMAX 
C  RESPECTIVELY. 

C 

C  IF  IRSET*1  IN  COMMON  /GS/.  THE  POINTERS  WILL  BE  RESTORED  ON  EXIT  TO 
C  THEIR  ORIGINAL  VALUES  WHEN  GAUSS4  WAS  ENTERED.  THIS  MEANS  THAT  THE 
C  SAME  GAUSSIAN  SEQUENCE  < EXCEPT  FOR  THE  WEIGHTS  A)  WILL  BE  GENERATED  ON 
C  THE  NEXT  CALL. 

C 

C  PRIOR  TO  CALLING  GAUSS4,  THE  RANDOM  NUMBERS  IN  ARRAYS  TR  AND  TI  MUST 
C  BE  INITIALIZED  WITH  A  CALL  TO  GAUSSI. 

C 

DIMENSION  A(U»  XR(1)»  XI(1) 

COMMON  /CB/  IRSET.  NTR.  NTI.  NMAX. KR, KI,  TR< 1253), TI (1353) 

IF(N.  GT.  NMAX  )  STOP 
DO  35  I-l.N 
XR(I)*A(I)*TR(KR) 

XI ( I )=A( I )*TI (KI ) 

KR-KR+1 
KI-KI+1 
35  CONTINUE 

IF( IRSET.  NE.  1)  GO  TO  40 

KR-KR-N 

KI-KI-N 

40  IF(KR.  GT.  NTR)  KR=KR-NTR 
IF(KI.  OT.  NTI)  KI-KI-NTI 
RETUBN 
END 


